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We present an update on the meson mass spectra calculation with the light-front quark model 
constrained by the variational principle for the QCD-motivated effective Hamiltonian. By smearing 
out the Dirac delta function in the hyperfine interaction, we were able to avoid the negative infinity 
that one encounters when using variational principle for the entire Hamiltonian. We obtained a 
better fitting for the mass spectra of ground state pseudoscalar and vector mesons from 7r to T, 
compared to the previous calculation that handled the hyperfine interaction as a perturbation rather 
than including it in the parameterization process. Our new parameters are further tested in decay 
constant calculation. We showed that by taking a larger harmonic oscillator basis in our trial wave 
function, the decay constants calculated from our model can be improved to agree reasonably well 
with the experimental data. 



I. INTRODUCTION 

Effective degrees of freedom to describe a strongly in- 
teracting system of hadrons have been one of the key 
issues in understanding the non-perturbative nature of 
QCD in the low energy regime. Within an impressive 
array of effective theories available nowadays, the con- 
stituent quark model has been quite useful in providing 
a good physical picture of hadrons just like the atomic 
model for the system of atoms. Absorbing the compli- 
cated sea quark pairs and gluonic interactions into the 
masses of constituent quarks and effective interactions in 
the Hamiltonian, one can make the problem much sim- 
pler yet still keep some key features of the underlying 
QCD to provide useful predictions [1]. The commonly 
used valence constituent quark models truncate the Fock 
space to the lowest sector and simplify the whole sys- 
tem to a bound state of dressed valence quark/antiquark 
constituents with QCD inspired effective potentials. Al- 
though a clear link between the QCD and this type of 
constituent quark models is still pending, the justifica- 
tion of effective potentials are usually provided by the 
flux tube configurations generated by the gluon fields as 
well as the one-gluon-exchange calculation in QCD [2, 3]. 

In addition, relativistic effects play an important role 
in describing the low-lying hadrons made of u, d, and 
s quarks and antiquarks. Thus, a proper way of deal- 
ing with the relativistic aspect of the hadron system is 
also quite essential. The formulation of light-front dy- 
namics(LFD) [4], or the front form dynamics introduced 
by Dirac [5] , provides a natural framework to include the 
relativistic effects which are crucial in describing the low- 
lying mesons. The simple vacuum except the zero-mode 
in LFD is a remarkable feature and facilitates an effective 
description of hadrons in a Fock space. In this formal- 
ism, each hadron is characterized by a set of Fock state 
wave functions, i.e. the probability amplitudes for find- 
ing different combinations of bare quarks and gluons in 
the hadron at a given LF time t = t + z/c. Moreover, the 
LF wave function is independent of all reference frames 



related by the front-form boosts because the longitudinal 
boost operator as well as the LF transverse boost opera- 
tors are all kinematical. This is clearly an advantageous 
feature unique to LFD, which makes the calculation of 
observables such as form factors much more effective than 
any other forms of relativistic dynamics. 

Taking advantage of the LFD, we have developed [6, 7] 
a light-front quark model (LFQM) based on a simple 
QCD-motivated effective Hamiltonian for a description 
of mesons. The key idea of our LFQM was to treat the 
radial wave function as a trial function for the variational 
principle to the QCD-motivated effective Hamiltonian 
with the well-known linear plus Coulomb interaction. 
Both the masses and the hadronic properties of ground 
state pseudoscalar and vector mesons in our LFQM were 
fairly well reproduced by taking just a lS'-state harmonic 
oscillator (HO) wave function as a trial function. Com- 
puting the meson mass spectra [6-8], however, we have 
treated the hyperfine interaction as a perturbation rather 
than including it in the variation procedure to avoid the 
negative infinity from the Dirac delta function contained 
in the hyperfine interaction. This treatment, however, 
left a room for any further development in handling the 
hyperfine interaction in our LFQM since the justification 
of using the perturbation was too technical to be backed 
up by the first principle, although our numerical results 
were quite comparable with the available experimental 
data not only for the meson mass spectra but also for 
the meson wave function related observables such as the 
decay constants, form factors, etc. 

The main purpose of this work is to update our pre- 
vious model [6-8] by including the hyperfine interaction 
also in our variational method to get the optimal model 
parameters, and examine if it improves our numerical re- 
sults compared to the ones obtained by the perturbative 
treatment of the hyperfine interaction. To achieve this 
goal, we smeared out the Dirac delta function by a Gaus- 
sian distribution in order to resolve the infinity problem 
when variational principle is applied. Moreover, we im- 
proved our trial wave function by taking a larger HO ba- 
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sis to analyze this effect of improving trial wave functions 
on our numerical results. To examine our updated model 
prediction, we calculate both the meson mass spectra and 
the meson decay constants. 

The paper is organized as follows: In Sec. II, we de- 
scribe our QCD-motivated effective Hamiltonian with the 
smeared-out hyperfine interaction. Using two different 
radial wave functions, i.e. the ground state HO wave 
function and the mixture of the two lowest order HO 
states, as trial wave functions of the variational principle, 
we find the analytic formulae of the mass eigenvalues for 
the ground state pseudoscalar and vector mesons. The 
optimum values of model parameters are also presented 
in this section. In Sec. Ill, we present our numerical re- 
sults of the mass spectra obtained from both trial wave 
functions and compare them with the experimental data 
as well as our previous calculations [6-8]. To test our 
trial wave functions with the parameters obtained from 
the variational principle, we also calculate the meson de- 
cay constants and compare them with the experimental 
data as well as other available theoretical predictions. 
We show an improvement of our numerical results by 
taking a larger HO basis in the trial wave function over 
the case of taking just the ground state HO wave func- 
tion as the trial wave function. Summary and conclusion 
follow in Sec. IV. The detailed procedure of fixing our 
parameters through variational principle is presented in 
the Appendix A. 



where (S q Sq) — 1/4 (—3/4) for the vector (pseudoscalar) 



meson and V 2 y cou i = 
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5 3 (r). The LF wave function 



of the ground state mesons is given by 

^(a^ku^) =^ J f_(a; i ,k_ L i)0(a; i ,k_ L i), (3) 



where <p is the radial wave function and Tt J \ J l_ is the 
interaction-independent spin-orbit wave function. The 
wave function is represented by the Lorentz invariant in- 
ternal variables Xi = p\ / P + , kj_i = pj_j — XiP ± and hc- 
licity X i: where = (P + , P~, P ± ) = (P° + P 3 , (M 2 + 
Pj_)/P + , P_i_) is the momentum of the meson with mass 
M and pf is the momenta of constituent quarks. 

The spin-orbit wave functions for pseudoscalar and 
vector mesons are given by [8, 9] 



oo 

A Aq 



n 



V2M 



(4) 



where ^{ J z ) is the polarization vector of the vector me- 
son and Mo = \/Mq — (m q — nig) 2 . The boost invariant 
meson mass squared Mq obtained from the free energies 
of the constituents is given by 
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II. MODEL DESCRIPTION 



In our LFQM for mesons, we approximate the sys- 
tem as effectively dressed valence quarks governed by the 
following QCD-motivated effective Hamiltonian in which 
motion of the quarks inside a meson is relativistic [6-8] : 



H = 



k 2 



l| + k 2 + V, 



(1) 



where k = (kj_,fc z ) is the relativistic three-momentum 
of the constituent quarks in the rest frame of the meson. 
To describe the interaction between quark and antiquark, 
we use the linear confining potential V^onf plus the one- 
gluon exchange potential Vogc. For S- wave pseudoscalar 
and vector mesons, the one-gluon exchange potential re- 
duces to the coulomb potential V^ ou i plus the hyperfine 
interaction Vh yp - Thus, we have 

V = Konf + v OKC 



(2) 




The spin-orbit wave functions satisfy the relation 
\ %\ J \'R\ J \ = 1 for both pseudoscalar and vector 
mesons. So the computation on this part of the wave- 
function is rather trivial in our approach. We will thus 
mainly focus on the radial wave functions in the follow- 
ing. 

To use a variational principle, we take our trial wave 
function as an expansion of the true wave function in the 
HO basis. We use the same trial wave function <\> for both 
pseudoscalar and vector mesons, but we try two different 
forms: one simply takes the form of the l^-state HO 
wave function given by 



iis(xi,k_ Lj ) 
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and the other one is expanded with the two lowest order 
HO wave functions 



V/ ^is(^i,kxi) + \/l - f cf>2s(xi,k±i), (6b) 
where / gives the fraction of 15* wave function and 
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For convenience, we will call Eq. (6a) and Eq. (6b) our 
wave function 4>a and (/>b, respectively. In these equa- 
tions, /3 is the variational parameter and k z is the longi- 
tudinal momentum defined by k z = (x — 1 /2) Mq + (to? — 
TOg)/2Mo. Thus the Jacobian of the variable transforma- 
tion (x,k±) —> k = (kj_, k z ) is given by 
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The normalization factor in <j) n s is obtained from the 
following normalization of the wave function: 

J' dxJ^\<f> nS (x t ,k ±l )\ 2 = l. (9) 

Next, we evaluate the expectation value of the Hamil- 
tonian in Eq. (1) with 0a(b)j i- e - (4 > A(B) \H\4>a(b)) which 
depends on the variational parameter j3. According to 
the variational principle, we can set the upper limit of 
the ground state's energy by calculating the expectation 
value of the system's Hamiltonian with a trial wave func- 
tion. In our previous calculations [6-8], which we call "CJ 
model" , we first evaluate the expectation value of the cen- 
tral Hamiltonian T + V con [ + V cou \ with the trial function 
<f>A, where T is the kinetic energy part of the Hamilto- 
nian. Once the model parameters are fixed by minimiz- 
ing the expectation value {4>a\{T + V con{ + Ko U i)|</>a), 
then the mass eigenvalue of each meson is obtained as 
M q q = (<f>A\H\(f)A) ■ The hyperfine interaction Vh yp in CJ 
model, which contains a Dirac delta function, was treated 



as perturbation to the Hamiltonian and was left out in 
the variational process that optimizes the model param- 
eters. The main reason for doing this was to avoid the 
negative infinity generated by the delta function as was 
pointed out in [10]. Specifically, (</>a| Vh yp |</>A) for pseu- 
doscalar mesons decreases faster than other terms that 
increase as /3 increases. So the expectation value of the 
Hamiltonian is unbounded from below. 

But now we want to include the hyperfine interaction 
in our parameterization process. To avoid the negative 
infinity, we use a Gaussian smearing function to weaken 
the singularity of S 3 (r) in hyperfine interaction, viz. [10, 
11], 

* 3 ( r ) ^72 e ^ V - ( 10 ) 

Once the delta function is smeared out as in Eq. (10), a 
true minimum for the mass occurs at a finite value of /3. 
The analytic formulae of mass eigenvalues for our modi- 
fied Hamiltonian with the smeared-out hyperfine interac- 
tion, i.e. M g ^ B ' = (<I>a(b)\H\4>a(b)) i are found as follows 

« yfRp V^P ^ V 2/3V 

v v i— q,q v / 

f 8/3 32(S g ■ S q )P 3 cr 3 

" [V^ 9m g m 5 0r(^ 2 + (T 2 ) 3/2 J' 

(11) 
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where i^i is the modified Bessel function of the second kind and U(a,b,z) is Tricomi's (confluent hypergeomet- 
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ric) function. 

We then apply the variational principle, i.e. 
dM^ B ^ /dj3 = 0, to find the optimal model parameters 
in order to get a best fit for the mass spectra of ground 
state pseudoscalar and vector mesons (a more detailed 
description of this procedure can be found in Appendix 
A). 

Our optimized potential parameters are obtained as 



{a 



-0.5575 GeV, b = 0.18 GeV z 



for and {a = -0.6664 GeV, b 



and a s = 0.5174} 
: 0.18 GeV 2 ,a s = 
0.5348} for M^, respectively. We should note that the 
two sets of potential parameters are quite comparable 
with the ones suggested by Scora and Isgur [12], where 
they obtained a = -0.81 GeV, b = 0.18 GeV 2 , and a s = 
0.3 ~ 0.6. For a comparison, the coupling constant we 
found in our previous model [6-8] was a s — 0.31. For the 
best fit of the mass spectra, the mixing factor / for <pB 
is found to be 0.7, i.e. 70% of the 15 state contribution 

to <fi B . 

Since we included the hypcrfine interaction with smear- 
ing function entirely in our variational process, we now 
obtained the two different sets of f3 values, one for pseu- 
doscalar and the other for vector mesons, respectively. 
The optimal Gaussian parameters j3 qq - for pseudoscalar 
and vector mesons are listed in Table I and II, respec- 
tively. Our optimal constituent quark masses and the 
smearing parameter a are also listed in Table III. 

Our modified model with the smeared-out hyperfine 
interaction improves the mass spectrum fitting, which is 
presented in the next section. This may suggest that 
when using constituent quark models, the contact in- 
teractions has to be smeared out. In fact, we think 
this smeared-out interaction seems to be more consis- 
tent with the physical picture for a system of finite-sized 
constituent quarks. 

For practical application of our model, we also compute 
the decay constants for the ground state pseudoscalar and 
vector mesons. The decay constants are defined by 



(0|fy"75<z|P) = i/pP", 
'\yfq\V(P,h)) = fvM v e>>(h), 



(13) 



for pseudoscalar and vector mesons, respectively. The 
experimental values of the pion and rho meson decay 
constants are f n rj 131 MeV from tt — > pv and /„ ps 
220 MeV from p -> e+er . 

Using the plus component (/i = +) of the currents, one 
can easily calculate the decay constants. The explicit for- 
mulae of pseudoscalar and vector meson decay constants 
are given by [8, 9] 
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where A = (1 — x)m q + xuiq and M.q = Mq + m q + rriq. 



FIG. 1. (color online). Fit of the ground state meson masses 
[MeV] with the parameters given in Table I, II and III for <pA 
and (j)B , compared with the fit from our previous calculations 
using CJ model [8] as well as the experimental values. The 
(•7T, p) masses are our input data. 



Here only the L z = S z = component of the wave func- 
tion contributes. Note that the vector meson decay con- 
stant fv is extracted from the longitudinal (h — 0) po- 
larization. We perform the decay constant calculations 
for both trial wave functions (/>a and using the corre- 
sponding set of parameters fixed for each trial wavefunc- 
tion, respectively. 



III. RESULTS AND DISCUSSION 

We show in Fig. 1 our prediction of the meson mass 
spectra obtained from the variational principle to the 
modified Hamiltonian with the smeared-out hyperfine 
interaction using two different trial functions 4>a (blue 
lines) and (purple lines) and compare them with the 
experimental data (green lines). We also include the re- 
sults (black lines) obtained from the CJ model with the 
linear confining potential [8]. We should note that the 
masses of 7r and p mesons are used as inputs in our 
calculation. As one can see, the single IS state HO 
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TABLE I. The Gaussian parameter j3 [GeV] for ground state pseudoscalar mesons obtained by the variational principle, 
q = u and d. 



Model 


Pqq 


Pqs 


Pqc 


/3cs 


/?cc 








/?bb 


0A 


0.6376 


0.5513 


0.5810 


0.5994 


0.7916 


0.6686 


0.7132 


1.0577 


1.6455 




0.4520 


0.3799 


0.3960 


0.4078 


0.5286 


0.4461 


0.4757 


0.6891 


1.0549 



TABLE II. The Gaussian parameter /3 [GeV] for ground state vector mesons obtained by the variational principle, q = u and d. 



Model 


Pqq 


Pqs 


Pqc 


Pes 


Pee 


Pqb 


Pbs 


p bc 


Pbh 


0A 


0.3480 


0.3952 


0.5283 


0.5727 


0.7849 


0.6436 


0.7010 


1.0554 


1.6450 


0B 


0.2416 


0.2742 


0.3579 


0.3892 


0.5233 


0.4278 


0.4671 


0.6871 


1.0544 



TABLE III. Constituent quark masses [GeV] and the smear- 
ing parameter a [GeV] obtained by the variational principle 
for the Hamiltonian with a smeared-out hyperfine interaction. 
Here q = u and d. 



Model 




m s 


m c 




a 


0A 


0.220 


0.432 


1.77 


5.2 


0.405 


0B 


0.221 


0.456 


1.77 


5.2 


0.423 



wave function a already generates a good enough fit- 
ting for the spectrum, and a more complicated trial wave 
function B does not change the IS results too much. 
In fact, the x 2 value for this modified model is 0.014 
[0.018] for 0a [0b], which is more than half reduced from 
\ 2 = 0.039 for the CJ model [8]. Except for the mass 
of K, our predictions for the masses of 15-state pseu- 
doscalar and vector mesons are within 4% error. Es- 
pecially, our modified Hamiltonian clearly improves the 
predictions of heavy-light and heavy quarkonia systems 
such as {r] c ,J/ij),B c ,r}b,T) compared to the CJ model 
adopting the contact hyperfine interaction. Although the 
experimental data for B* is not yet available, our predic- 
tions of B*, i.e. 6343 (6325) MeV for A(B ), are quite 
comparable with other quark model predictions such as 
6345.8 MeV [13] and 6340 MeV [11]. 

In Table IV, we list our predictions for the decay con- 
stants of light mesons (w,K,p,K*) obtained by using 
0A(B) and compare them with CJ model [14] and the 
experimental data [15, 16]. As one can see, our updated 
model calculation including the hyperfine interaction in 
the variation procedure doesn't seem to improve the re- 
sults of CJ model. In particular, the trial wave function 
0A generates decay constants that are quite high for light 
mesons (it, p, K, K *) indicating that just lS'-state HO 
wave function alone cannot be a good trial wave func- 
tion for the entire Hamiltonian including the smeared 
hyperfine interaction. However, using the improved trial 
wave function 0b, we can see a dramatic decrease in the 
numerical results consistent with the variational princi- 
ple. Indeed, the results from 0b are much closer to the 



experimental data than those from 0a- Especially for n, 
the decay constant changes from 155 MeV to 139 MeV, 
which is much closer to the experimental value. Although 
the CJ model yields the experimental value of / w much 
better than the updated model results from 0b, an over- 
all improvement due to the change of trial wave functions 
from 0a to B seems quite clear. Since the experimental 
values are very well known for light mesons, this improve- 
ment is very encouraging. 

In Table V, we list our predictions for the charmed 
meson decay constants (f D , f D * , f Ds , f D * , f Vc , fj/y) to- 
gether with CJ Model [17], lattice QCD [18, 19], QCD 
sum rules [20], relativistic Bethe-Salpeter (BS) model 
[21], relativized quark model [22], and other relativistic 
quark model (RQM) [23] predictions as well as the avail- 
able experimental data [15, 16, 24, 25]. We extract the 
experimental value (/j/*) C xp = (416 ± 6) MeV from the 
data r exp (J/* -> e+e") = 5.55 ± 0.14 ± 0.02 keV and 
the formula 

4-7T rv 2 

V(V e+e") = YMy^cv, (15) 

where cy = 4/9 for V — J/^. We should note that our 
results of the ratios JdJId = 1.13 [1.14] and f^/fjfq, = 
0.88[0.91] obtained from 0a [0b] are quite comparable 
with the available experimental data, /d s //d = 1.25 ± 
0.06 [15, 16] and f^Jfj/a = 0.81 ±0.19 [24, 25], respec- 
tively. Our results of the ratios /d*//d* = 1.15[1.17] 
obtained from 0a [0b] are also in good agreement with 
the BS model prediction, f D *JfD* = 1-10 ± 0.06 [21] 

We list our results for the bottomed mesons 
(/s> /s* , /s a , /b* , /»)„ , /t ) in Table VI, and compare with 
CJ Model [17], lattice QCD [18, 26, 27], QCD sum rules 
[20, 28], BS model [21], relativized quark model [22], 
and RQM [23] predictions as well as the available ex- 
perimental data [24, 29]. Note that we extract the ex- 
perimental value (/y)cxp = (715 ± 5) MeV from the 
data r cxp (T -> e+e ) = 1.340 ± 0.018 keV [24] and 
Eq. (15) with cy = 1/9 for V = T. Our results for 
the ratios f B Jf B = 1.17(1.18], f Bi /f B * = 1.19[1.20] ob- 
tained from 0a [0b] are in good agreement with the re- 
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cent lattice results, 1.20(3)(1) [26] and 1.22 (tg) [27] for 
f B Jf B and 1.17(4)1] [18] for /b ; // b .. Our results for 
the ratio f Vb /fr = 0.92 [0.95] are to be compared with 
the f Vb / fr ~ 1 in HQSS limit [30]. One can also see that 
for heavy charmed and bottomed mesons, the trial wave 
function 0b produces better results when compared with 
the experimental data as well as the the lattice results. 

In Table VII, we present our model predictions for the 
decay constants of fs c and fs* , and compare them with 
other model calculations [8, 31-36]. Both trial wave func- 
tions a and </>b give results that are comparable with 
other model calculations. 



tonian. Our calculation seems to suggest that our modifi- 
cation of the Hamiltonian by smearing out the hyperfine 
interaction contributes to an improvement in the right 
direction. 

Since our modified Hamiltonian together with the im- 
proved trial wave function provides very good results as 
discussed in this work, it may be also desirable to in- 
vestigate further and see if we can improve our previous 
calculations of other wave function related observables 
such as form factors. We shall explore them in our fu- 
ture work. 



IV. SUMMARY AND CONCLUSION 
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In this work, we updated our LFQM by smearing out 
the Dirac delta function in the hyperfine interaction and 
including the smeared hyperfine interaction in our calcu- 
lation based on the variational principle rather than us- 
ing the perturbation method to handle the delta function 
in the contact hyperfine interaction. Using the two trial 
wave functions, i.e. the IS" state HO wave function 0a 
and the mixed wave function 0b of IS* and 25* HO states, 
we calculated both the mass spectra of the ground state 
pseudoscalar and vector mesons and the decay constants 
of the corresponding mesons. We found that our predic- 
tions of the meson mass spectra are in good agreement 
with the data both for a and 0b, which seems to sug- 
gest that for the energy eigenvalues, a simple trial wave 
function like 0a may suffice. However, the mixed wave 
function 0b of IS and 2S states turns out to be much 
better in the calculation of decay constants as the wave 
function a generates results that are noticeably higher. 
According to the variational principle, the results from 
0b are obviously closer to the true results compare to 
the results from 0a- This improvement is especially ev- 
ident for the light mesons (n, K, p, K*) where the decay 
constants are very well known. 

The fact that by including one more higher order HO 
basis in the wave function can improve the results for the 
decay constants significantly is quite encouraging. Since 
the more realistic wave function for our modified Hamil- 
tonian in Eq. (1) can in principle be written in a complete 
HO basis, i.e. J2"=T c n4>nSi the improvement obtained 
by 0b — XiLi c n<t>nS is consistant with our expectation. 
Nonetheless, we should point out that, strictly speaking, 
we are not comparing trial wave functions for exactly 
the same Hamiltonian here, since we not only varied the 
parameters in our wave functions, but also adjusted the 
parameters in the Hamiltonian, such as the masses of 
quarks and the smearing parameter for best fitting of 
the mass spectra, as can be seen in Table III. This makes 
the Hamiltonian we used for 0a and 0b slightly differ- 
ent. However, the parameters in the Hamiltonian are 
kept within a limited range and thus one can still antici- 
pate that a more sophisticated trial wave function should 
produce better predictions from the same form of Hamil- 
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Appendix A: Fixation of the model parameters 
using variational principle 

In our model, we assumed SU(2) symmetry and have 
the following parameters that need to be fixed: con- 
stituent quark masses (m u r,n, m s , m C) mb), potential pa- 
rameters (a, 6, a s ), gaussian parameter /?, and the smear- 
ing parameter a. For our second trial wave function 0b, 
we also have the mixing factor / that we have to adjust. 
Notice that the j3 values here are not only different for 
different quark combinations, but also different for pseu- 
doscalar and vector mesons of the same quark combina- 
tion. The reason for this is that the hyperfine interaction 
we included in our parameterization process gives differ- 
ent contributions to the masses of pseudoscalar and vec- 
tor mesons and thus induces different parameterizations 
under variational principle. 

We now illustrate our procedure for fixing these param- 
eters. The variational principle gives us one constraint: 



d{&\H\y) dM q 



Oft 



dp 



o. 



(Al) 



We can use this equation to rewrite the coupling con- 
stant a s in terms of other parameters and plug it back 
into Eq. (11) [Eq. (12)] and thus eliminate a s . The string 
constant b is fixed to be 0.18 GeV, a well known value 
from other quark model analysis [11, 12, 37]. We will 
leave the quark masses and smearing parameter a (and 
the mixing factor / for 0b) as externally adjustable vari- 
ables. We picked a set of values for (m u /,n , m s , m c , nib, a) 
when using 0a and (m u u), m s, ni c , nib, &, f) when using 
0b, and proceed with the following procedure to solve for 
the rest of parameters. We then vary these values to do 
the procedure again until we find a set of quark masses 
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TABLE IV. Decay Constants for light mesons (in unit of MeV) obtained from our updated LFQM . 



fir fp Ik fK' 

Model <j>A 155 234 190 261 

Model <j) B 139 211 176 234 

CJ Model [14] 130 246 161 256 

Exp. [15, 16] 130.4 ±0.2 221 ± 1 156.1 ±0.8 217 ±7 



TABLE V. Charmed meson decay constants (in unit of MeV) obtained from our updated LFQM. 





fo 


fo- 


fo. 


fo: 






Model A 


244 


279 


276 


322 


406 


460 


Model B 


218 


241 


249 


282 


354 


390 


CJ Model [17] 


197 


239 


232 


273 


326 


360 


Lattice [18] 


211 ± 3± 17 


245 ± 20i^ 


231 ± 12±l 


272 ± 16:^ 






QCD [19] 


201 ±3± 17 




249 ± 3 ± 16 








Sum-rules [20] 


204 ± 20 




235 ± 24 








BS [21] 


230 ± 25 


340 ± 23 


248 ± 27 


375 ± 24 


292 ± 25 


459 ± 28 


QM [22] 


240 ± 20 




290 ± 20 








RQM [23] 


234 


310 


268 


315 






Exp. 


206.7 ±8.9 [15, 16] 




257.5 ±6.1 [15, 16] 




335 ± 75 [25] 


416 ± 6 [24] 



TABLE VI. Bottomed meson decay constants (in unit of MeV) obtained from our updated LFQM. 

!b !b- !b 3 !b* fr lb /t 

Model A 229 243 267 288 805 871 

Model B 195 202 229 242 654 692 

CJ Model [17] 171 185 205 220 507 529 

Lattice [18] 179 ± 18±| 4 196 ± 24±f 9 204 ± let^ 1 229 ± 20^ 

QCD1 [26] 216 ± 22 259 ± 32 

QCD2 [27] 189 ± 27 230 ± 30 

Sum-rulesl [28] 210 ± 19 244 ± 21 

Sum-rules2 [20] 203 ± 23 236 ± 30 

BS [21] 196 ± 29 238 ±18 216 ± 32 272 ± 20 498 ±20 

QM [22] 155 ±15 210 ± 20 

RQM [23] 189 219 218 251 

Exp. 229±|f±|^ [29] 715 ± 5 [24] 



TABLE VII. Bottom-charmed meson decay constants(in unit of MeV) obtained from our updated LFQM. 





4>A 




CJ Model [8] 


[31] 


[32] 


[33] 


[34] 


[35] 


[36] 


/Be 


488 


406 


349 


360 


433 


500 


460 ± 60 


517 


410 ± 40 




530 


432 


369 




503 


500 


460 ± 60 


517 





that give best fit for the meson mass spectra. In addi- /3^, are the gaussian parameters for pseudoscalar (p) 

tion, we have 3 more parameters (a,/3^,P^) to be fixed and vector (v) mesons, respectively. 

for mesons of a certain quark combination fog), where TT • f A , n t 

^ W1 " Using the masses ot 7r and p as our input values lor 



MA in Eq. (11) [Mf 9 in Eq. (12)], and the condition 
that our coupling constants a s are the same for all these 
ground state pseudoscalar and vector mesons, we can fix 
the three model parameters (a, f3 qq , (3 qq ) for q = u or d 
from the following three equations: 

M w (!3 p qq ,a) = 0.140, (A2a) 
M^-, a) = 0.77, (A2b) 
a a (J% 9 ,o) = a.(% q ,a). (A2c) 

Solving these equations not only gives us the remaining 
parameters a, 0t q , and/3",, but also the coupling constant 



a s which we assumed to be the same for all the mesons 
we consider here. We can then solve for the f3 values of all 
the other mesons using the known a s value, by equating 
the a s expressions for different mesons that we got from 
Eq. (Al). We thus fixed all the parameters for all the 
mesons we consider here. 

Through our trial and error type of analysis, we 
found m q = 0.220 GeV,m s = 0.432 GeV,m c = 
1.77 GeV,m b = 5.2 GeV,cr = 0.405 GeV gives best fit 
of the meson mass spectrum when using trial wave func- 
tion 4>ki while m q = 0.221 GeV, m s — 0.456 GeV, m c 
1 .77 GeV, nit, = 5.2 GeV,a = 0.423 GeV, / = 0.7 gives 
best fit when using trial wave function </>b- For these 
values, our obtained (3 values are listed in Table I and II. 
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